Assessment of 5-year outcomes of life satisfaction in survivors after rehabilitation programs: a multicenter clinical trial

Using a rehabilitation program for the survivors of acute respiratory distress syndrome (ARDS) could be one of the important and fundamental steps to improve the pulmonary function and health-related quality of life (HRQoL) of patients. This study was carried out to evaluate the effect of two rehabilitation techniques (Family-Based Empowerment Model (FECM)/Continuing Care Model (CCM), or both of them) on pulmonary function, and HRQoL in ARDS survivors. From December 2009 to June 2016, ARDS survivors from mixed medical-surgical ICUs at four academic teaching hospitals in Tehran, Iran, were randomly assigned to one of three intervention groups (A, B, or C) or a control group (D). Pre- and post-interventions, pulmonary functions and HRQoL status of patients in all groups were collected 48 times via clinical measurements and various questionnaires during 5 years of follow-up. Significantly improvement was seen in the intervention groups compared to the control group, and the greatest benefit was observed in patients who received mixed of FCEM and CCM rehabilitation techniques. Co-administration of FCEM and CCM can improve pulmonary function as well as the life satisfaction of ARDS survivors. As a result, the execution of the empowerment model by nurses is recommended for ARDS survivors and the participation of their families at the same time. Trial registration: NCT02787720 (ClinicalTrial.gov, 24/05/2016).

, and the familycentered empowerment model (FECM) 14 . CCM model focuses on the effective, interactive, and balanced roles of healthcare providers, the patient, and the support structure so as to effectively influence the patients' attitude toward the disease, improve treatment compliance and participation. This model was associated with improved HRQoL in renal transplant recipients 15 , patients following coronary artery bypass graft surgery 16 , and the patients with diabetes mellitus 17 .
FCEM was developed to improve the care and outcomes of the patients with chronic diseases and has previously been evaluated and validated in several chronic disease states 14 . Family engagement in the rehabilitation process may potentially have beneficial effects, including the improvements in mental and physical function, and reduced burden of disease 18,19 . The primary aim of model is to empower the patient/family unit to promote health quality. The model has four stages: (a) determining perceived threat (group discussion method); (b) self-efficacy (problem-solving method); (c) improving self-esteem (educational participation method) and (d) process and outcome evaluations 14 .
Little data is available comparing the effectiveness of pulmonary rehabilitation models in the ARDS survivors. We investigated the impact of pulmonary rehabilitation programs on ARDS survivors using either FCEM and CCM models or a hybrid FCEM-CCM model as compared to routine care on HRQoL.

Results
Participants of the study. From December 2009 to June 2016, 283 out of 367 ARDS survivors in the mixed medical-surgical ICU at four academic teaching hospitals in Tehran, Iran, met the inclusion criteria. Two-hundred eighty-three eligible patients were included and evenly randomized into four groups. A total of 143 patients due to reasons such as permanent ventilator dependence (n = 21), change of code status from "full code" (n = 27), revoked consent (n = 23), death (n = 24), and incomplete data (n = 48), were excluded from the final analysis. Therefore, 35 patients remained in each group, and 140 patients entered the final analysis. Figure 1 shows the flowchart of participants in the trial study.
Demographic and clinical characteristics at baseline and over 5-years. 140 subjects were included in the analysis. The mean age of the participants was 62 ± 8 years, with more than half of the patients was female (n = 82, 58.1%). Demographic and clinical characteristics at the baseline are presented in Table 1 for all groups of study. Patient's demographics were similar among groups, including age (P = 0.165), sex (P = 0.525), BMI (P = 0.590), pre-ARDS BMI ≥ 30 (P = 0.24), ARDS risk factors (P = 0.85), family size (P = 0.12), marriage status (P = 0.156), urban residence (P = 0.63), and working full-time or part-time (P > 0.05). In addition, no significant difference was observed among the groups of study in terms of clinical characteristics. The only exception was in renal replacement therapy, which was significantly different among the groups (P = 0.043). Table 2 are shown the results of repetition a 5-year summary of demographic and clinical variables, including organ dysfunction, pulmonary dysfunction, coexisting illness, smoking status, and ability to return to work   Table 2). All 3 groups as intervention groups (groups A, B and C) were merged and compared to the control group (group D) in both models as the additional analysis. Return to work differed significantly between intervention and control groups, with both GEE and Panel analysis showing a higher chance of returning to work for group A than groups B, C, and D. The likelihood of returning to original work was greatest for Group A, followed by groups B, C, and D. The results held true for the groups A, B, and C combined vs. group D as controls.
The odds ratio (OR) for GEE and Panel models indicated a decrease in coexisting illness for patients in intervention groups vs. controls, with group A < B < C < D, respectively. No significant association between existing organ dysfunction and intervention group was noted (P > 0.05), but the chance of existing pulmonary dysfunction was statistically lower for group A vs. controls (group D) based on GEE (OR: 0.38, P = 0.021) and Panel analysis (OR: 0.007, P = 0.027). No significant differences were noted for the groups B and C (P > 0.05). Although smoking was different in the fourth and fifth years of follow-up, only the GEE model indicated a significant difference between group A and controls (OR: 0.430, P = 0.04) (Supplementary file, Fig. 2).
Health-related quality of life findings. Table 3 presented the results for HRQoL variable scores, includ-ing QoL, perceived stress, state anxiety, trait anxiety, Barthel index (BI), Kessler psychological distress scale (K10), six-minute walk test (6MWT), and a free walking index (WI) test. QoL improved via the study across all groups, with significantly greater improvements noted in the intervention groups (A, B, and C) compared to controls (61.3 ± 4.7 vs. 24 ± 6, P < 0.0001). The greatest benefit was observed in Group A that the mean score of QoL was significantly increased from 21 ± 2 to 86 ± 3 (P < 0.0001). GEE and Panel analysis models (Supplementary file, Table 3) indicated positive effects of the intervention on increasing BI via 5-year follow-up for intervention groups vs. control group (17 ± 1.8 vs. 10 ± 1.8, P < 0.0001). During the study, K10 scores decreased across all groups with significantly greater improvements noted in the intervention groups (A, B, and C) compared to controls (24 ± 6 vs. 33.5 ± 4, P < 0.0001), and the greatest benefit was observed in group A (23 ± 6 vs. 26 ± 5, P < 0.0001). However, state anxiety increased across intervention groups (A, B and C) as compared to group D, with the greatest increase belonging to group A (45 ± 7 vs. 80 ± 0.2, P < 0.0001). Trait anxiety remained stable   www.nature.com/scientificreports/ throughout the study in groups B, C, and D (P > 0.05), while decreasing in group A (54.5 ± 4 vs. 54 ± 3, P = 0.042).
6MWT index increased for all groups during the 5-year follow-up. Significant improvement was observed for aggregate intervention groups (A, B, and C) compared with controls (518 ± 213 vs. 418 ± 260, P = 0.032); however, the greatest individual benefit was observed in group A (71 ± 5 vs. 533 ± 221, P < 0.0001). The WI also increased for all groups throughout follow-up. When comparing aggregate intervention groups (A, B and C) to controls, significant improvement was noted in the intervention cohort (P < 0.05). However, when individually comparing groups A, B, and C to group D, only group A was significantly better than controls (23,202 ± 8040 vs. 18,871 ± 8421, P = 0.025). Pre-and post-intervention, as well as during a 5-year follow-up, repeated measurements of life satisfaction items were evaluated by soft outcome via two SF-36 questionnaire items (physical and mental components score) and two hard outcomes (K10) and (BI) in pre-and post-intervention and during the 5-year follow-up (Fig. 2). Besides, to determine the time of better rate of BI and K10 index over 5-year follow-up, we used a multiple Cox survival analysis (Supplementary file, Table 1). The results for BI indicated better scores for intervention groups, compared to controls as the reference group. In group A, the hazard ratio (HR) of the score of BI upper than 14 was (HR: 87.65, P < 0.0001) compared to controls (Supplementary file, Fig. 2a). Furthermore, the probability of a higher score of BI (> 14) increased for those who reported only one coexisting illness through the follow-up (HR: 2.045, P = 0.017). Besides, the HR of the lower K10 index was (HR: 2.179, P = 0.048), and (HR: 2.409, P = 0.009), in the first and second intervention groups, respectively. While it was not significant for third groups (HR: 0.615, P = 0.264), compared to controls (Supplementary file, Fig. 2b). Besides, BMI (P = 0.139) and Returned to work first year (P = 0.170) indicated no prognosis on KPDS. Table 4 shows the outcomes of pulmonary function testing on ARDS survivors during a 5-year period. Throughout the follow-up period, all groups' total lung capacity (TLC) rose. GEE and Panel analysis indicated that TLC improvement amongst combined intervention groups (A, B, and C) was significantly improved in comparison with controls (91 ± 7 vs. 83 ± 9, P < 0.0001). Individual improvements were noted for group A > B > C, respectively. The same improvement patterns were also observed for forced expiratory volume in 1 s (FEV 1 ), forced vital capacity (FVC), and the diffusing capacity of lungs for carbon monoxide (DLCO). The ratio of FEV1 to FVC (FEV1/FVC), was another significant increased pulmonary function test which increased in all intervention groups, compared to controls (80 ± 10 vs. 73 ± 6, P < 0.0001), according to GEE and Panel analysis (Supplementary file, Table 4), but for total intervention groups, it was not significant (P = 0.109). Finally, RV test was decreased via the 5-year follow-up for intervention groups vs. controls and total intervention groups according to both models (P < 0.0001). SEM findings. SEM was performed to identify direct and indirect factors influencing outcomes in the patients experiencing. The normed χ 2 was 1.57, indicating excellent fit. Moreover, the RMSEA = 0.064, IFI = 0.99, and CFI = 0.99 indicating that the model, including these factors performed better to describe the data. Therefore, the modified model fits the data (Supplementary file, Table 5).

Pulmonary function findings.
Among many factors included in the SEM to evaluate the impact of clinical factors on ARDS patient outcomes, independent factors, including TLC, DLCO, FVC, FEV 1 , FEV 1 /FCV, WI, and (6MWT) were significantly increased by pooled intervention group (groups 1-3). While the effect of this variable decreased RV and K10. DLCO directly correlated with K10 (β = 0.208, P < 0.05), and K10 displayed a significant positive correlation with   Fig. 3).

Discussion
This study was conducted to determine the effect of FCEM and CCM rehabilitation techniques on pulmonary functions and HRQoL of ARDS survivors. According to the results, overall life satisfaction, including QoL, mental health (perceived stress and psychological distress), and physical health (physical state, activities of daily living, and exercise capacity), improved during the study in all groups, and significantly improvement was seen in the intervention groups compared to the control group, and the greatest benefit was observed in patients that used the mixture of FCEM and CCM rehabilitation techniques (group A). However, state anxiety increased among the intervention groups compared to the control group, and the highest increase was related to the first group that received FCEM and CCM programs, and trait anxiety was stable during the study in groups B and C as intervention groups and group D as the control group, while in the first group (group A) a decrease was www.nature.com/scientificreports/ observed. The stable situation of trait anxiety and the increasing rate of state anxiety can be a reflection of the patient's sense of empowerment in managing his or her own health. Besides, our interventions do not include forced treatment or interventions that can change the nature of patients' anxiety. Therefore, this level remained stable 14 . In terms of pulmonary function, measured factors, including TLC, FEV 1 , FVC, FEV 1 /FVC and DLCO increased and RV decreased in all groups throughout the 5-years follow-up. Analysis indicated that the pulmonary factors significantly improved among intervention groups compared to the controls. A few studies reported the recovery of pulmonary functions, exercise capacity, and QoL in patients who survived ARDS20-22. A study by Hsieh et al. 11 , showed improvement in pulmonary function and exercise capacity in ARDS survivors of influenza A (H1N1) at 3 months after hospital discharge due to exercise, simple strength training, and respiratory exercise as SPR. Although, they reported the improvement of the QoL of these patients at 6 months after hospital discharge, even though there was no further improvement of their pulmonary functions and exercise capacity. Evidence suggests that lung volume shows a strong tendency to return to normal 3 to 6 months after the acute phase 20,21 . However, 6% to 43% of patients develop an obstructive pattern, and 15% to 58% develop a restrictive pattern in the first year of follow-up 22 . In addition, previous studies showed that 1 year after discharge, the 6-min walking distance has increased compared to the findings of immediate postdischarge period 9,23 .
To date, no studies are examining the effect of FCEM and CCM methods and comparing them to standard pulmonary rehabilitation on ARDS survivors. However, few studies were performed on the effects of these methods on the patients with other chronic illness [15][16][17][18][19] . The results of these studies about using an FCEM on caregivers and family members of patients with chronic illness, have been shown some improvement in the quality of life and mental health, which was consistent with our findings and reported the implementation of an FCEM intervention can reduce the burden of disease, and improve the mental and physical health of patients [24][25][26] . According to them, FCEM reduces the burden of care through its three axes: motivation, psychology (self-confidence, selfcontrol, and self-efficacy), and problem-solving capacity. In terms of effect CCM, several studies showed that nurses could apply the CCM as an effective method to reduce risk factors and improve the lifestyle of patients with chronic disease30. However, a study conducted by Mohammadi et al. 27 , to determine the effect of using CCM at home on the QoL of patients with myocardial infarction (MI) showed that a rehabilitation program consisting of training sessions on MI and its complications, diet food, and medicine, risk factors, etc. at home did not have a significant effect on different dimensions of QoL of patients and there was no significant difference among the groups. These discrepancies could be related to sample size, patient follow-up period, and type of chronic disease.
In general, the present study results showed that in all groups, in the post-intervention phase, pulmonary function, exercise capacity (6MWD and WI), and HRQoL improved. These improvements were significantly higher in the intervention groups than in the control group. Among the intervention groups, the most improvement occurred in group "A" patients who used FCEM and CCM rehabilitation methods. The findings support the hypothesis that the synergistic effect of FCEM and CCM methods can significantly improve the patient's physical and mental function, reduce the burden of disease and improve the pulmonary function of ARDS survivors. CCM model's major purpose is to provide ongoing follow-up treatment in order to manage the illness and any potential consequences. Over the course of 12 weeks, four steps of identification, sensitization, control, and assessment are completed; the first two take 3 weeks while the other two take 9 weeks. The FCEM model's primary purpose, however, is to empower the family system and the patient as a whole in order to enhance their health. Perceived danger (severity and sensitivity), self-efficacy, self-esteem (confidence), and appraisal are the phases. As a result, the CCM model's basic and crucial concept of "follow-up and continuity of care" and the FCEM model's key concept of "the inseparable whole of the family and patient system" were merged; the synergistic effect of the two very important concepts mentioned had a dramatic effect on the first group's results when compared to the other two intervention groups (each model separately). Besides, encouraging the patient,  www.nature.com/scientificreports/ understanding the patient more, reinforcing positive and feeling responsive via the simultaneous participants of family or friends (FCEM), and being in contact with a nurse (CCM), which creates a greater sense of support and satisfaction. This is reflected in the results of physical and mental health scores. To our knowledge, this is the first research that evaluates and compares the synergistic effect of FCEM and CCM methods on pulmonary functions and HRQoL variables in ARDS survivors with standard pulmonary rehabilitation. Besides, this multi-center study was well-designed with a suitable sample size and 5-year followup. Nevertheless, this study had some limitations, which related to its nature. It was a randomized, controlled, and blinded prospective study, and 5 years of follow-up of patients was very difficult and in terms of the loss of follow-up has led to a decrease in sample size.
Finally, co-administration of FCEM and CCM may enhance pulmonary function as well as ARDS survivors' life satisfaction. As a consequence, ARDS survivors and their families are encouraged to participate in the empowerment paradigm, which should be carried out by nurses. Besides, further studies regarding the synergistic effect of FCEM and CCM models on the care performance of caregivers, as well as the pulmonary functions and QoL of ARDS survivors due to side effects, are suggested. On June 1, 2016, the experiment was registered with Clinicaltrials.gov (NCT02787720). On release from the ICU, the patient provided written permission. In circumstances where the patient lacked decision-making ability, surrogate consent from the patient's legal guardian or healthcare proxy was allowed. Furthermore, the research followed the guidelines of the Helsinki Declarations 29 .

Methods
Eligibility criteria of participants. From December 2009 to June 2016, ICU patients with ARDS were screened for eligibility. The diagnosis of ARDS was based on the Berlin definition 1 , which was as follows; patients with a BMI < 40 suffering from ARDS with PaO 2 /FiO 2 < 300 mmHg during mechanically ventilated (MV) with an expected duration of controlled MV of more than 24 h and ability to tolerate PEEP titration (up to 21 or 15 cmH 2 O). Besides, the patients were eligible to enroll in this study if they met the following criteria: (a) age ≥ 18 years, (b) able and willing to provide informed consent, (c) willingness of a family member or friend of patients to participate in the study, (d) has basic health literacy and can fill out questionnaire, and (e) full code status. Patients aged less than 18 years and more than 85 years, pregnant women, patients at end-stage medical condition, patients with a history of pulmonary rehabilitation, and patients with a history of neurological or psychiatric disorders were excluded from the study.
Sample size, randomization and blinding. According to power estimates, each group would need 32 patient/family units to obtain a 95 percent confidence level and 90% power. Convenient sampling was used to enlist participants. In terms of sample loss, the sample size in each of the intervention and control groups was eventually determined to be 35. Eligible patients admitted to ICU were enrolled within 24 h and randomly assigned to one of three intervention groups (A, B, or C) or a control group (D). Group (A) who has received both FCEM and CCM programs via trainer. Group (B) received only FCEM via family member/friend as a trainer. Group (C) who received only CCM via researcher as trainer, and group (D) who used routine care as a control group. Randomization was achieved with a computer-generated random block design, which was drawn up by an expert statistician who had no clinical involvement in the trial before the beginning of study. Block randomization was accomplished using Random Allocation Software© (Informer Technologies, Inc., Madrid, Spain) by computer-generated random. Randomly allocated numbers were placed into sequential containers (i.e. ABCD, containers 1-4; BADC, containers 5-8, etc.) which were kept in a secure location until allocation consignment. The difference among groups was not disclosed, patients consented knowing that they were undergoing pulmonary rehabilitation (PR), but without knowing the details. The assignment was made through confidential communication between the patient's and a third party not involved in the recruitment process. Hence, patients and their family members, as well as data analyzers, were blinded to the assignment group and the differences among the groups.
Intervention. The intervention package had three phases: pre-intervention, intervention, and post-intervention phases (Supplementary file, Fig. 1, and additional explanation about methods).
Role of the designee and rehabilitation plan. Following informed permission, the selected family member or friend (hence referred to as designee) remained with the patient as a 'unit' throughout the trial. During stages 3 and 4, the designee attended the patient's educational sessions, with stage 2 being optional based on the designee's preferences. All patients had similar inpatient rehabilitation programs. For the patients in three intervention groups, outpatient rehab included daily exercise for 0-2 h/day. Exercise occurred between 8:00 and 10:00, and types included walking, jogging, bicycle, swimming, or other exercises according to patient preference or resource availability and confirmed by the multidisciplinary medical team.  Statistical analysis. The sample size in each group was calculated to be 32, with a confidence level of 95 percent and a test power of 80 percent, assuming that the measure of the impact of FCEM and CCM on parental burden of care (middle effect size) is at least d = 10. In terms of sample loss, the sample size in each of the intervention and control groups was eventually determined to be 35. In order to control dropouts in this trial, 140 participants were recruited. Statistical analysis was performed using IBM ® SPSS ® Statistics 21.0 (IBM Corp., Armonk, NY) and IBM ® SPSS ® AMOS™ 21.0 (IBM Corp., Armonk, NY). Discrete variables are expressed as counts and percentages. Power calculations determined that 32 patient/family units were needed in each group to achieve a 95% confidence level and a 90% power. Continuous variables are expressed as means and standard deviations. Analysis of variance (ANOVA) and Chi-Square analyses was used to compare numeric and discrete variables. Generalized Estimation Equations (GEE) and Panel Analysis were performed on longitudinal data, and the results were expressed as odds ratios (ORs). Multiple Cox proportional hazards regression was used to assess the effect of treatment groups on Barthel time and KPDS time. Structural Equation Modeling (SEM) (more information about SEM are available in Supplementary file) was performed to examine the effects of clinical factors on outcomes of acute respiratory distress syndrome, and the model was evaluated using the root mean square error of approximation (RMSEA), the normed fit index (NFI), and the goodness-of-fit index (GFI). Statistical significance was defined as p-value < 0.05.

Ethical approval. The protocol study was reviewed and approved by the Ethics Committees of Hamadan
University of Medical Sciences (IR.UMSHA.REC.1400.604). The trial was registered with Clinicaltrials.gov (NCT02787720) on 24/05/2016. In addition, the study was conducted according to the Helsinki Declarations guideline.

Data availability
All data collected and analyzed during the current study are available from the corresponding author on reasonable request.